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Abstract 

In this paper an approach is proposed to introduce Hamiltonian formalism to dissi- 
pative mechanical systems. This approach is based upon below viewpoints proposed 
in this paper: under a certain identical initial condition a dissipative mechanical sys- 
tem shares a common phase flow curve with a conservative system; the Hamiltonian 
of the conservative system is the value of the total energy of the dissipative system 
under the initial condition; whether the dissipative mechanical system has an ana- 
lytical solution or not, the conservative system may always exist, which is to say the 
existence of the Hamiltonian and Hamilton's equation. Hence this approach is to 
substitute an infinite number of conservative systems for a dissipative mechanical 
system corresponding to varied initial conditions. The major means to demonstrate 
these viewpoints is that by the Newton-Laplace principle the nonconservative force 
can be reasonably assumed as a function of a component of generalized coordinates 
qi along. The advantage of this approach is such that there is no need to change the 
definition of canonical momenta and the motion is identical as that of the original 
system; this approach is suitable to all dissipative mechanical systems regardless of 
the existence of analytical solution. 
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1 Introduction 



Since Hamilton originated Hamilton equations of motion and Hamiltonian 
formalism, it has been stated in most classical textbooks that the Hamiltonian 
formalism focuses on solving conservative problems. 



In 1960s, iHori and Brouwerl (Il96ll ) utilized the classical Hamiltonian formal- 
ism and a perturbation theory to solve a non-conservative problem. They did 
not attempt to derive t he Ha miltonian formalism for non-conservative prob- 
lems. In 1990s iTveterl (119981 ) started an attempt and thereby obtained his 
general Hamilton equation: 



p- = -(—) +f(—) 

\ d Qi) Q Pt \ d Qi) Q Pt 

where {Q, P} denote the canonical variables and the position vector r depends 
on the canonical variable set {Q, P} and t, i.e. r(Q, P,t), K is the transformed 
Hamiltonian and the subscript in the partial derivative expressions indicate the 
functional dependencies of K and r. If the variable set {Q, P} is transformed 
into the variable set {q,p}, where p is mechanical momentum, the position 
vector r depends on q and does not depend explicitly on t, i.e., r(q), Eq. (11.11) 
could be reduced to 



Pi 



\dqij qp \dq t 




X2) 



where F(dr/dqi) q denotes a generalized force in direction i and H is the 
Hamiltonian. Yet both Eq. ( 11. II) and Eq. ( 11.21) are not truly Hamiltonian 
formalisms, because there is not a conservative Hamiltonian quantity (first 
integral). In 1990s other similar attempts were launched, but the resulting 
formalisms were not truly Hamiltonian. 



In 1940s IP. Morse and Feshbachl (119531 ) gave an example of an artificial Hamil- 
tonian for a damped oscillator based on a "mirror-image" trick, incorporat- 
ing a second oscillator with negative friction. The resulting Hamiltonian is 
un-physical: it is unbounded from below and under time reversal the oscilla- 
tor is transformed into its "mirror-image". By this arbitrary tri ck dissipa- 



tive systems can be handled as though they were conservative. IVujanovic 
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(119891) invest i gated dissipative systems from the vision of variational meth- 
ods. iTarasovl (120011 ) suggested a generalization of canonical q uantizati o n that 
maps a dynamical operator to a dynamical superoperator. ITarasovl (120011 ) 
claimed that this approach allows defining con sistent quan t izatio n procedure 



for non-Hamiltonian and dissipative systems. iR.M.Kiehnl (119741 ) considered 
that dissipation effects may be included by considering the dissipative sys- 
tems for which the closed integral of action is a parameter-dependent, confor- 
mal in variant of the motion. He applied this idea to hydrodynamics. iRajeev 
( 120071 ) considered that a large class of dissipative systems can be brought to 
a canonical form by introducing complex coordinates in phase space and a 
complex-valued Hamiltonian. The papers all took into consideration, the en- 
ergy drained from the dissipative system. In general they used an example 
of a damped oscillator to demonstrate their approach. But the equations de- 
scribing the oscillator always were one-di mensional or dec oupled to a group of 
one-dimensional equations. For instance, IVujanovid (119891 ) offered an example 
: for the linear damped oscillator mx + fix + cx = 0, the following conservation 
law can be obtained: ((x 2 /2) + (uj 2 x 2 /2) + kxx)e 2kt = const, where u 2 = c/m 
and 2k = fi/m. We can see that e 2kt is a decay coefficient of amplitude. We 
think, if and only if the linear damped oscillator is one dimensional or can 
be decoupled to a group of one-dimensional equations, the decay coefficient 
can be presented as e 2kt . In addition we consider, that if one were to apply to 
these approaches to dissipative systems, one must first change the definition 
of canonical momenta and canonical coordinates. 



Based on Tveter's Eq. (11.21) . in this paper we attempt to apply Hamilton for- 
malism to mechanical systems with damping forces. The approach is to substi- 
tute a group of conservative systems for a dissipative mechanical system. We 
have found that one can take the value of the total energy of the dissipative 
system as the Hamiltonian of a conservative system, the motion of which is 
the same as that of the dissipative system with a certain initial condition. 



2 The derivation 



2.1 Obtaining the Hamiltonian Quantity of the substituting conservative sys- 
tem 



Under general circumstances, the force F is a damping force that depends on 
the variable set qi, ■ ■ ■ , q n , q%, ■ ■ ■ , q n . We denote by Fj the components of the 
generalized force F. 



Fifai, ■ ■ ■ , Qn, Qi, ■ ■ ■ , Qn) = F — . (2.1) 
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Thus we can reformulate the Eq. fll.2p as follows: 



(9H\ 

Pi = ~ J + /•;■(</!• ■ • • • </„.</)• • ' ••'/») 



(2.2) 



'//< 



Suppose the new Hamiltonian Quantity is H and we do not change the defi- 
nition of the canonical momenta, thus we may write a Hamilton's equation of 
a new conservative system: 



Pi = ~ 



. d( ii , 



'IP 



(2.3) 



where the initial condition is same as that of Eq. fl2.2p . i.e. p Q = Po,Qo — %■ 
Under the initial condition described as above, we assume p = p, q = q. By 
comparing Eq. (12.21) with Eq. fl2.3p . we have 



In classical mechanics the Hamiltonian H is mechanical energy: 



(2.4) 



(2.5) 



where 7 denotes a phase flow curve presented by Eq. (12.21) . c\ is a constant that 
depends on the initial condition above. The Einstein summation convention 
has been used in this paper. 

Hence an attempt is made to find H through line integral along the curve 7: 



(©.*- / 1 



'dH\ 
fa) 



F i (q 1 ,---,q n ,q 1 ,---,q n ) 



up 



dq~i 
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m iPi= [( a JL) dp,. (2.6) 



v <>i>i , , , 

7 N 'IP 7 'IP 

In like manner with Eq. (l2.5p . we have 

H=l{ d S) ^+f($) #, + c 2 , (2.7) 



j \ / qp j \ •■ / qp 



where c 2 is a constant which depends on the initial condition. On the path 7 
we have 

dqi = dqi,dpi = dp (2.8) 



Substituting Eq. (l2.5l)(l2.6p(l2.8p into Eq. fl2.7p . it can be represented as: 

H = H - J Fi{q u ■ - ■ , q n , q u ■ ■ ■ , q n )dqi + c. (2.9) 
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where c = c 2 — ci. According to the physical meaning of Hamiltonian, c\ and c 2 
and c are added into Eq. (12. 5p (12.61) (12.91) such that the integral constant vanishes 
in the Hamiltonian quantity. It is well known that Eq. (12.21) denotes a phase 
flow curve, thus according t o the Newt o n-Lap lace principle of determinacy 



described in book written by I V.I. Arnold! (119781 ). we can assume 



4i = <ii(t), 

where the flow curve satisfies the initial condition. Hence we divide time do- 
main into a group of sufficient small domain and consider a small domain, 
such that we can reasonably assume F{ as 

Fi(qi(t(qi)), ■ ■ • , g n (*(*)), ?i (*(?.)), • • ' , ?n(*(?i))) = ^(ft), 



where Ti is a function of qi alone. Indeed T is a conservative force, the value 
of that is equal to the value of F% at 7. Thus we have 



J Fid qi = J F^dqi = - W^). (2.10) 
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According to the physical meaning of Hamiltonian, c is added to Eq. (12.91) 
such that the integral constant vanishes in Hamiltonian quantity. Hence c = 
-Wi(qf). Substituting Eq. ([2TTUj) and c = -Wi(qf) into Eq.Q), we have 

H = H-W l (q i ) (2.11) 



where — Wi(qi) denotes the negative work done by the damping force F. In 
Eq. (12.1 II) H and H are both functions of g^ and Wi(qi) a function of gj. If we 
substitute and pi for g^ and p i; and do not change the form of functions H 
and H, Wi(qi), we would then obtain the Hamiltonian of Eq. (12.31) 

H{q h Pi) = H{q u Pi)Wi{qi) (2.12) 



Then we must show that the Hamiltonian presented by Eq. (l2.12l) satisfies the 
Eq. (12.41) under the same initial condition, i.e. that the motion of Eq. (l2.2p is 
equivalent to that of Eq. (l2.3p under the same initial condition. 

Substituting Eq.flmj into the left side of Eq.flMD, we have 



dH(qi,pi) _ dH(q h pi) dWj(qj 



dH(qi,pi) = dH(qi,pi) dWjjqj) 
dpi dpi dpi 

It must be emphasized that although g« and p. t are considered as distinct 
variables in Hamilton's mechanics, we can consider g^ and q\ as dependent 
variables in the process of constructing of H. At the trajectory 7 we have 



dpi 



0, (2.14) 



where Fi(qi) is the value of damping force Fi on the phase flow curve 7. 
Substituting Eq. (l2.14p into Eq. (12. 131) and substituting g^ and Pi for g^ and pi, 
Eq. (12. 131) can be reduced to Eg. (12.41) under the initial condition. 

Therefore we can say that the motion of Eq. (l2.3p is equivalent to that of 
Eq. (12.21) at the initial condition; the dissipative system presented by Eq. (12.21) 
can be replaced by the conservative system presented by Eq. (l2.2p at the initial 
condition; and H is the Hamiltonian of a substituting conservative system. 
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2. 2 Example 



Take an oscillator with a spring and damping as an example, the governing 
equation of which is as below: 

q + cq + kq = 0, (2.15) 



where the stiffness coefficient k is a constant, and the damping coefficient c 
also is a constant. The solution of Eq. (12.15j) is 



q = Ae~^ cos ( {C2 - Ak) \ + B 



2 

,2 



(2.16) 



q = _^(c 2 - Ak)e- Ct ' 2 sin f ( ° 4fc) % + B 



-A- °-e'^ cos (^plt + B^ (2.17) 

where constants A and B are determined by an initial condition. We assume 
t belongs to a sufficient small domain, so we let 

t = t(q) (2.18) 



As the section above by substituting Eq. (12.18j ) into Eq. fl2.17| ) one can represent 
the damping force function of q alone: 

cq = cq(q), (2.19) 



which can be rewritten as: 

cq = k(t)q. (2.20) 

Hence one obtains a new conservative system as below: 

q + (k + k{t))q = 0. (2.21) 

This is a nonlinear oscillator, the stiffness of which varies with the time t. In 
order to make it easy to distinguish Eq. fl2.21D from Eq. fl2.15p . we denote the 
position variable in Eq. fl2.2ip as q. If and only if the Equation above has the 
same initial condition as Eq.f l2.15l) . Eq. (12.211) shares the same solution with 
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Dissipative Mechanical Systems with varied Initial Conditions 




Conservative System 1 



Conservative System 2 



Conservative System n 



Fig. 1. 



Dissipative Mechanical System and Conservative Systems 

Eq. (12.151) . The Hamiltonian quantity of Eq. (12.211) is the value of the total 
energy of the Eq. (l2.15p under the initial condition: 



H = ~q 2 + htq 2 + / rqiXq. 



[2.22) 



According to Eq. ( 12. 19112. 20p . the Hamiltonian can be represented as: 
H=\q 2 + \kq 2 + J k(t)qdq = -q + -kf + J k(t)qdq 

7 7 



(2.23) 



It must be emphasized that !^k{t)qdq is a function of q alone, and k(t) is 
determined by the initial condition of Eq. (12. 151) . As such the Hamiltonian 
quantity H is determined by the initial condition of Eq. (l2.15l) . 



2.3 Discussion 



Based on the above, we can outline the relation between a dissipative me- 
chanical system and conservative systems with Fig([T]). As shown in Fig([T|), 
the relation between a dissipative mechanical system and conservative sys- 
tems can be stated as below: 
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The initial condition of the original dissipative system determines the phase 
flow cure 7; the value of the total energy of the original system H depends on 
7; one can take H as the Hamiltonian of a substituting conservative system; 
the substituting system shares only one common phase flow curve 7 with the 
original system. 

Although the substituting system shares a common phase flow curve with the 
original system, under other initial conditions they exhibit different phase flow 
curves. Therefore the phase flow of the substituting systems differs from that 
of the original s ystem . It is well known that the Liouville theorem as explained 



by IV.IArnoldl (119781 ). can be formulated as: 



Theorem 2.1 The phase flow of Hamilton's Equations preserves volume: for 
any region D we have 

volume of g l D = volume of D 



where g l is the one-parameter group of transformations of phase space 
^:(p(0),g(0))^: (p(t),q(t)) 



According to the theorem above, the phase flow of the original dissipative 
system Eq. (12.21) certainly does not preserve its phase volume, but the phase 
flow of the substituting conservative Eq. (12.31) does. 

The Hamiltonians of new conservative systems in general are not analytically 
integrable, unless the original mechanical system is integrable. The reason is 
that the work done by damping force depends on the phase flow curve. If 
the system is integrable, then the phase flow curve can be explicitly written 
out, the system has an analytical solution, and therefore the work done by 
damping force can be explicitly integrated. Therefore the Hamiltonian H can 
be explicitly represented. Most systems do not have an analytical solution. 
Despite this, the Hamilton quantity and coordinates and m omentum must 



satisfy the Eq. (l2.3p under a certain initial condition. Why had IF. Klein! (119281 ) 
written: "Physicists can make use of these theories only very little, an engineers 
nothing at all"? The answer: when one is seeking an analytical solution to a 
classical mechanics problem by utilizing Hamiltonian formalism, in fact one 
must inevitably convert the problem back to Newtonian formalism. This means 
that an explicit form of Hamiltonian quantity is not necessary for classical 
mechanics. What is important is the relation between q,p and Hamiltonian 
quantity embodied in the Hamilton's Equation. 
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3 Conclusion 



We can conclude that a dissipative mechanical system has a property: under 
a certain common initial condition the dissipative system shares a common 
phase flow curve with a conservative system; the Hamiltonian of the conser- 
vative system is the value of the total energy of the dissipative system. Based 
on this point it can be concluded, that a dissipative problem can be refor- 
mulated as an infinite number of non-dissipative problems, corresponding to 
every phase flow curve of the dissipative problem. One can avoid having to 
change the definition of the canonical momentum in the Hamilton formalism, 
because under a certain initial condition the motion of one of the group of 
conservative systems is same as the original dissipative system. 

Because conservative systems may own some other fine characteristics (e.g. 
capable of preserving symplectic structure), it is feasible to substitute a group 
of conservative systems for the original dissipative mechanical system, and 
thereby introduce the Hamiltonian formalism to dissipative mechanical sys- 
tems. 
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